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ABSTRACT 

High-resolution iV-body simulations of hierarchical cosmologies have shown that 
the density and velocity dispersion profiles of dark-matter haloes display well-definite 
universal forms whose origin remains unknown. In the present paper, we calculate the 
internal structure of haloes expected to arise in any such cosmologies by simply tak- 
ing into account that halo growth proceeds through an alternate sequence of discrete 
major mergers and long periods of gentle accretion. Major mergers cause the violent 
relaxation of the system subject to the boundary conditions imposed by accreting lay- 
ers beginning to fall in at that moment. Accretion makes the system develop inside-out 
from the previous seed according to the spherical infall model. The predicted structure 
is in very good agreement with the results of numerical simulations, particularly for 
moderate and low mass haloes. We find strong indications that the slight departure 
observed in more massive systems is not due to a poorer theoretical prediction, but 
to the more marked effects of the limited resolution used in the simulations on the 
empirical profiles. This may have important consequences on the reported universality 
of halo structure. 
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1 INTRODUCTION 



Hoffman 200( ; Del Popolo et al. 200C ; Engineer, Kanekar & 



Padmanabhan 2000; Nusser 2001) 



The dominant dark component of matter in the universe 
appears to be clustered in bound haloes which form the 
skeleton of all astronomical objects of cosmological inter- 
est, from dwarf galaxies to rich galaxy clusters. All we know 
about these systems comes from their gravitational effects 
on the luminous matter they trap or on the light traveling 
across them. This is so little information that the following 
fundamental questions have prevailed for a long time. Does 
the internal structure of haloes depend on their mass? Does 
it depend on their past history? And on cosmology? How is 
this structure set? 



An important result of this research line is that haloes 
grow inside-out as new material is being incorporated in sec- 
ondary infall. That is, despite the continuous shell-crossing 
(and radial relaxation, see below) of the infalling layers with 
the previously relaxed body, the structure remains at any 
moment essentially unaltered within the instantaneous ra- 
dius. As shown below, this sole consideration should per- 
mit us to determine the halo density profile provided the 
infall rate of matter is known. Much progress has been 
achieved in the last twenty years in the modelling of halo 
growth. The extended Press - Schcchter (PS) model ( Press & 



k. n^tr nQ79; |n»tr ^Q7^ |n,.r.r, iq7-;| ) w^r^ tv.^ 



Schechter 1974 



first tola,ddress this problem by considering the simplified 
case of haloes forming through spherical infall, that is, the 
collapse of a density fluctuation of smooth, isotropically 
distributed, dissipationless matter in an otherwise homoge- 
neous expanding universe. Under the adiabatic invariant col- 
lapse approximation, they derived the density profile arising 
from some speciflc initial conditions. The effects of refining 
the derivation used and of adopting more and more realis- 
tic initial conditi ons have been subsequently addressed in 
a seri es of works (Fillmore fc Goldreich 1984 ISertschingei 



Bower 1991 



Bond et al. 1991 



Lacey fc Cole 



199; ) makes quite accurate prediction s, indeed, on the rate 



19851; iHoffman fc Shaham 198E|; iRyden fc Gunn 1987|; H 



den 1988t [Zaroubi fc Hoffman 1993| ; pLokas 200C| ; [Lokasfe 



at which haloes increase their mass (Lacey & Cole 1994) 
in hierarchical cosmologies like the one describing the real 
universe. Unfortunately, haloes develop in such cosmologies 
through continuous mergers rather than through smooth 
spherical infall. 

The effects of mergers depend on the relative mass of the 
progenitors. Major mergers bring the whole system out of 
equilibrium making it move in the phase space around some 
attractor. During this process, particles experience random 
accelerations owing to the rapidly varying collecti ve poten- 
tial well which causes the relaxation of the system ( Lynden 
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Bell 19 67) in a similar although more dramatic way than 



in the case of spherically infalling layers. This violent relax- 



ation erases any imprint of the previous history of the sys- 
tem, the new equilibrium state reached being characterised 
by a normal distribution of particle velocities similar to that 
yielded by the common two-body relaxation, although inde- 
pendent of particle mass. Were haloes isolated and, hence, 
unperturbed after any such dramatic event they would end 
up as infinite, spherical systems with a uniform, monomass, 
isotropic velocity dispersion. This is the reason why haloes 
are often r nodelled as isotropic, isoth ermal spheres (e.g.. 
King 1972; shapiro, Iliev & Raga 1999). But haloes are not 
isolated systems. During the violent relaxation process, they 
keep on collecting matter through minor mergers making the 
boundary conditions to vary in some unknown way. This 
severely limi ts the predict ive power of the violent relaxation 
theory (e.g., Saslaw 1987). 

To gain a deeper insight on the internal structure 
of haloes many authors have turned to numerical exper- 
iments. Using high-resolution cosmological A'^-body simu- 
lations Navarro, Frenk & White (1996, 1997) have found 
that the spherically averaged density profile of haloes of all 
masses is well fitted by the simple analytical expression 



p(r) 



r(rs + r) 



(1) 



in any hierarchical cosmology analysed. In equation (|l|), r 
is the radial distance to the halo centre, and and pc are 
the halo scale radius and characteristic density, respectively. 
The latter two parameters are related to each other and 
to the mass M of the halo through the condition that the 
virial radius R of the system encloses, by definition, an av- 
erage density equal to a fixed factor / times some density pu 
characterising the universe at that moment. The values of / 
found in the literature are in the range [178, 500], while pu 
is taken equal to pcrit, the critical density for closure, or p, 
the mean cosmic density. Navarro and collaborators adopted 
/ = 200 and pu = pcrit- 

On the other hand, the spherically averaged, locally 
isotropised, velocity dispersion profile S(r) is found to be 
well fitted by the solution of the Jeans equation for hydro- 
static equilibrium and negligible rotation. 



E^(r) 



dluE^ dlnp 



dlnr 



dlnr 



3GAf(r) 



(2) 



with p(r) the function given in equation (|l|) and M(r) the 
corresponding mass within r, by adopting the boundary con- 
dition of null pressure at infinity. Note that, although the 
hydrostatic condition is natural to hold, the boundary con- 
dition at infinity is not obvious owing to the limited extent 
of haloes. 

The se empirical results have been confirmed by other 



20017 



authors (Cole fc Lacey 1997; Huss, Jain & Steinmetz 1999b 



Bullock et al. 200l|), there being some controversy only at 



very small radii (fVIoore et al. 199^; |Jing fc Suto 2000| 



Fukushige fc Makino 2001) where they are the most affected 



by the spatial resolution of the simulations. In any event, the 
fundamental question about the origin of that empirical halo 
structure remains. It is not even clear whether this is a gen- 
eral result of gravitational collapse, including smooth spher- 
ical infall, or the specific consequence of repeated mergers 



(IHuss, Jain fc Steinmetz 19994 ^yer fc White 1998^ [Moore 
et al. 1999| ; [Lokas fc Hoffman 2000| 



In the present paper, we use a variant of the extended 
PS model distinguishing between minor and major mergers 
(§1^) to derive the structure of haloes predicted in hierarchi- 
cal cosmologies and compare it with that found in numerical 
simulations (§ The result of this comparison supports the 
idea that the internal structure of haloes is the natural im- 
print of their hierarchical growth (§0). 



2 MINOR AND MAJOR MERGERS 

The scaling of the previous density profile varies with cos- 
mology although the same trend is always found: the char- 
acteristic density pc decreases with increasing halo mass. 
Navarro, Frenk & White (1997) suggested that this univer- 



sal trend is due to the fact that, in hierarchical cosmologies, 
more massive haloes form later when the mean density of 
the universe is lower. Unfortunately, the idea that haloes fix 
their structural properties at formation is hard to assess be- 
cause the own concepts of halo formation and destruction 
are rather fuzzy in hierarchical cosmologies in which haloes 
grow through continuous mergers. 

To properly define these concepts it is con venient to 



distinguish between n ii nor and major mergers ( iManrique 
fc Salvador-Sole 1996|; [Kitayama fc Suto 19"9^; ^alvador- 



3ole, Solanes fc Manrique 1998~ | Percival, Miller fc Peacock 



200C; Cohn, Bagla fc White 2001). In the Modified Press- 



Schec hTer (MPS) model ( [Salvador-Sole, Solanes fc Manrique 



199S Rai^ j Gonzalez-Casado fc Salvador-Sole 1998 ; Raig 



Gonzalez-Casado fc Salvador-Sole 2001) mergers producing 
a fractional mass increase above a given threshold Am ~ 0.5 
are regarded as major, and minor otherwise. Note that a 
given merger may be seen as major or minor depending on 
the viewpoint of the partner halo which is considered. Since 
haloes are essentially unperturbed in minor mergers while 
they are completely rearranged in major ones we say that 
haloes survive in the former case and are destroyed in the 
latter. The destruction of a halo does not necessarily imply 
the formation of a new one: the largest partner partaking 
of the merger can see it as minor and, hence, survive to it 
being identified to the final system. This kind of mergers 
therefore correspond to the accretion by the most massive 
halo of the remaining partners. (The definitions of minor 
and major mergers given above imply that there is at most 
one surviving halo in any given merger.) Only those merg- 
ers in which all partners are destroyed or, equivalently, no 
merging halo survives, do mark the formation of new haloes. 

The MPS model yields analytical expressions for the 
mass accretion rate of haloes (see eq. ^ below), their de- 
struction and formation rates, as well as the distribution 
probability functions of progenitor masses and of formation 
and destruction times whose median values define, respec- 
tively, the typical mass of progenitors and the typical halo 
formation and destruction times. All these theoretical quan- 
tities arc in good agreement with the results of A^-body sim- 
ulations (Raig, Gonzalez-Casado fc Salvador-Sole 2001). 

In particular, the average rate at which mass is accreted 
onto a halo with M at t is, according to the definition of Am, 
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where 



nM(l+A„) 



{M' - M)rt{M, M',t)dM' , 



rt{M,M',t) dM' = 



/27^ d(5c d(j(M') 



a2(M') dt dM' 



a2(M) 



(3) 



,^-| -3/2 



X exp 



2cr2(M') 



(M') 



cr2(M) 



dM' 



(4) 



is the instantaneous transition rate at t from haloes with 
M to haloes between M' and M' + dM' due to mergers of 
whateve r amplitude provide d by the original extended PS 
model (Lacey & Cole 1993). In equation (0), 5c{t) is the 



linear extrapolation to the present time to of the critical 
overdensity of primordial fluctuations collapsing at t, and 
(j{M) = a{M,to) is the r.m.s. fluctuation of the density 
field at to smoothed over spheres of mass M; both Sc{t) and 
cr(M) depend on cosmology. 

The M{t) track followed during accretion by haloes with 
a given mass Mi at a given time ti is therefore the solution 
of the differential equation 



dM r,^/ N 1 

— ^r4M{t),t] 



(5) 



for the initial conditions M{ti) — Mi. Strictly speaking, this 
solution is the average track followed by those haloes. Real 
accretion tracks actually diffuse from it owing to the effects 
of individual random minor mergers, the scatter remaining 
non etheless quite limited along the typical lifetime of haloes 
(see |Raig, Gonzalez-Casado fc Salvador-Sole 2001 \ 



As shown by |Raig,^Son^Alez^Casado''ir'Salv^ 

(20010 the formation of haloes corresponds to rare binary 
mergers between similarly massive progenitors. Hence, they 
do cause the rearrangement of the corresponding global sys- 
tems, the resulting violent relaxation yielding a more or less 
spherical haloes with new density profiles independent of 
their past history. In contrast, minor mergers are frequent 
multiple events, particularly in the case of very small cap- 
tured haloes. The graininess of accreted matter can therefore 
be neglected in a first approximation, and the resulting con- 
figuration (i.e., a central more or less spherical, relaxed ob- 
ject surrounded by a rather smooth distribution of matter 
falling into it) can be well approximated by the spherical 
infall model. Consequently, haloes should grow inside-out 
during the accretion phase. 



3 THEORY VS. SIMULATIONS 

The point of view that haloes grow through long periods 
of gentle accretion after their violent formation in major 
mergers is clearly much more accurate than the two extreme 
points of view so far considered in order to derive analytical 
density profiles (see § |^), namely, the cases of pure spherical 
infall or one single isolated major merger. During the accre- 
tion phase, the inside-out growth condition holds and, since 
the average mass accretion rate is well-known (eq. ^), one 
can readily infer the structure developing around the initial 
seed. On the other hand, the structure of this latter, fixed 
at formation, should also be possible to determine since the 
boundary conditions met by the system during violent re- 



laxation are well-known, too: they are set by the accretion 
regime starting at that moment. 

3.1 Density profile 

Suppose a spherical halo formed at tf with mass Mf and 
radius R{. The typical M{t) track followed by the halo dur- 
ing the subsequent accretion phase is therefore the solu- 
tion of the differential equation (^ with initial condition 
M{t{) = Mf. According to the spherical infall model, the 
accreted mass is deposited, at any moment t, at the instanta- 
neous radius R{t) without altering the inner density profile. 
Consequently, we have 



M{t) - Ml = 



R(t) 



■iivr p{r) dr 



(6) 



Rl 



with p{r) the developing density profile at t. By differenti- 
ating this relation taking into account the abovementioned 
definition of the instantaneous radius 



3 M{t) 



1/3 



and equation ^ we are led to the expression 



dt 



(7) 



(8) 



r4M{t),t] 

for the density at r = R{t). Equations (^) and ^ therefore 
define, in the parametric form, the wanted density profile at 
any radius r ^ R{. 

But this is not enough to explain the shape of the den- 
sity profile drawn from simulations. Indeed, a significant 
fraction of it (of order 50 %) corresponds to the inner re- 
gion, r < Rf, fixed at formation (see the position of Rf in 
Fig. 2). To derive the density profile in that inner region we 
shall take into account that the structure emerging from vi- 
olent relaxation at tf must satisfy the boundary conditions 
imposed by accretion starting at that moment. This implies 
that the inner and outer profiles of any physical quantity 
must match up at R{. Some discontinuity would only ap- 
pear in case that these boundary conditions were incompat- 
ible with the inner structure. But, violent relaxation takes 
some time to proceed during which surrounding layers begin 
to be accreted. Hence, the final steady state, far from being 
incompatible with those boundary conditions, will perfectly 
adapt to them. 

Both M{t) and Pu{t) entering in equations (^) and (^ 
are known functions which can be differentiated to any arbi- 
trary order. Consequently, one can derive the values of p(r) 
and any order derivative of it at any radius in the outer 
region, in particular at the matching radius i?f. Then, by 
taking the Taylor series expansion of p(r) at Rf one can 
estimate that function at any radius r < 7?f to any arbi- 
trary accuracy. Although this proves the existence of one 
unique inner solution, a much practical way to derive it is 
the following one. As mentioned in § |^, all haloes having 
the same mass at a given time follow the same accretion 
track along the common time interval since their respective 
formation times tf, and since the functions M{t) and pu{t) 
coincide for all these haloes in that common time interval, 
the same is true for the respective outer profiles in the corre- 
sponding radial ranges. Moreover, the inner solutions of all 
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these haloes must coincide with the outer solutions of those 
having undergone the last major merger early enough since 
the profiles and derivatives of any order of these latter take 
identical values at the radii -Rf of the former and, there- 
fore, have identical Taylor series expansions at such radii. 
Strictly speaking, this conclusion refers to the average den- 
sity profile since individual profiles will show some scatter 
around it coming from the abovementioned scatter in real 
M{t) tracks. 

The previous reasoning proves from the mathematical 
point of view that all haloes with a given mass at a given 
time have the same (average) density profile regardless of 
their specific formation time, and that this unique density 
profile, which naturally appears after the past history of the 
system has been erased, coincides with the outer solution of 
haloes formed at an arbitrarily early epoch. From the phys- 
ical point of view the reason for this apparently surprising 
result is well understood. As well known, in the case of spher- 
ical infall there is only one mass distribution at any moment 
preserving the steady state of the inner system despite the 
shell-crossing of the infalling layers; this corresponds to the 
solution of the spherical infall model for a given accretion 
rate. In the case of violent relaxation subsequent to a ma- 
jor merger, the resulting steady structure must also adapt 
to the shell crossing of the infalling layers beginning to be 
accreted at that moment. This therefore necessarily leads to 
the same inner mass distribution as in the case of spherical 
infall at t hat moment. Note that this expla ins the result ob- 
tained by Huss, Jain & Steinmetz (1999a) that the density 
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profile of haloes closely resembles the NFW profile regard- 
less of whether they form in a hierarchical cosmology or via 
smooth spherical collapse. 

The exact values of / and pu used to delimit relaxed 
haloes in simulations or in the real universe are a mere con- 
vention and do not obviously affect their shape. However, 
from e quations (^ and (^) we see that the density pro- 



Figure 1. Theoretical density profiles for current haloes of masses 
IQ-^M*, IQ-^Mt, Mt, and 10 Af« (from bottom to top) in a flat, 
Ujn = 0.25, CDM model inferred using — Pcrit (dashed lines) 
and pu = p (full lines). In the plot we use the same delimitation 
of haloes as in Navarro et al. (1996, 1997). Af* is the mass scale 
at which primordial density fluctuations leave the linear regime, 
and p* stands for paitMt/M. 



ply reflects the real situation met in dealing with simulated 
(or observed) haloes: their actual steady region extends well 
beyond any popular delimiting radius. 

The correct behaviour of the extended PS model has 
only been checked in detail in the case of scale-free cosmolo- 



gies ( Lacey fc Cole 1994 Raig, Gonzalez-Casado & Salvador- 



file hero predicted depends on the definition of halo radius 



adopted. There is no contradiction however between these 
two facts. That definition is crucial, in general, for the ability 
of the PS model to correctly describe the halo growth pro- 
cess. For instance, the Press-Schechter mass function only 
fits the empirical mass function provided some specific defi- 
nition of halo radius as this fixes the mass of objects. Like- 
wise, we cannot pretend to recover the actual density profile 
of haloes from the PS formalism for any arbitrary definition 
of halo radius; this must be chosen so to obtain the best 
agreement between the predictions of the PS model and the 
results of simulations. Note however that such an "accret- 
ing radius" must not necessarily coincide with the arbitrary 
"delimiting radius" used to compare the resulting theoret- 
ical profile with the empirical one obtained by any specific 
author. In this latter case, we will operate with the theo- 
retical profile previously derived just as if it were the real 
profile of a simulated (or observed) halo: we shall apply the 
same delimiting radius which may or may not coincide with 
the accreting radius previously used. This can be readily 
done in most cases without any further consideration. Only 
in the case that the value of / used in the derivation of the 
profile is larger than the one used to delimit it in the final 
comparison, or that pcrit is used in the derivation while p 
is used in the comparison some extension of the theoretical 
profile towards the future will be required. But this sim- 



Sole 2001 for the MPS version). It is therefore necessary to 
establish the best value of pu to be used in case that pcrit 
and p do not coincide. In Figure 1, we plot the theoreti- 
cal profiles inferred from pu = pcrit and pu = p in a flat, 
Sim = 0.25 CDM cosmology. As can be seen, pu = p leads to 
a decreasing density profile with monotonous varying slope 
while pu = Pcrit leads to a density profile which tends to level 
off at large radii. (The radii where the two solutions devi- 
ate from each other correspond to cosmic times at which A 
becomes non- negligible.) Since such a behaviour is not ob- 
served in simulated (or observed) haloes, we conclude that 
the best value of pu to use is p. Concerning the value of / 
we have checked that the theoretical density profile is always 
very insensitive to it, any value in the range 178-500 yield- 
ing fully indistinguishable results at the resolution of Figure 
1. Hereafter, we adopt / = 200. 

In Figure 2, we compare the average density profiles 
predicted for current haloes of four different masses in three 
distinct cosmolog ies with the corresponding empi rical pro- 
files obtained by Navarro, Frenk & White (1997). As can 



be seen, there is good agreement between theory and sim- 
ulations, particularly for moderate and low mass haloes. A 
similar result is found in any cosmology analysed. We want 
to stress that, once we have fixed the value of the effective 
threshold Am for major mergers entering in the MPS model, 
the predicted density profiles of haloes of any mass at any 
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Figure 2. Predicted (full lines) and em pirical ("dotted lines') density pro files for current haloes with the same masses as in Fig. 1 in the 
three following cosmologies analysed by Navarro, Frenk fc White (1997 ): a flat, A = 0, CDM model (SCDM), a flat, Qm = 0.25, CDM 
model (ACDM), and a flat, A = 0, model with power- law spectrum of density fluctuations of index n = — 1 (FPL). Arrows mark the 
frontier between the inner regions fixed at the typical formation time and outer ones developed during the subsequent accretion phase. 
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Figure 3. Same density profiles as in Fig. 2 drawn using a value of the threshold for major mergers. Am, equal to 0.5 (full lines) 
compared to the profiles resulting from Am = 0.3 (dot-dashed lines) and Am = 0.7 (dashed lines). 



time are completely determined both in shape and normal- 
ization. T hus, there is no freedom to play with. As mcn - 
tioned in Raig, Gonzalez-Casado & Salvador-Sole (2001) 



the comparison between theory and simulations in terms of 
the description of halo growth alone does not allow us to 
constrain the value of Am. The only restriction is given a 



The observed trend for the agreement between theory 
and simulations to degrade towards large halo masses cannot 
be due to a bad modelling of the halo growth by means of the 
MPS model. As previously mentioned, the correct behaviour 
of the mass accretion rate has been checke d against numer 



ical simulations in scale-free c osmologies (R,aig, Gonzalez 



priori b y the expected dynamical effects of mergers of differ- 



Casado & Salvador-Sole 2001) such as the one plotted in 



ent amplitude which points to the plausible range [0.3, 0.7]. 
In Figure 3 we show the effects of adopting any of these 
two extreme values of Am. By comparing these results with 
those obtained from the most reasonable central value of 0.5 
shown in Figure 2 we conclude that the value of Am is, on 
the contrary, well-constrained around that central value by 
the internal structure of haloes. 



the right panel of Figure 2 where the departure is the most 
apparent. To understand the possible origin of this effect 
we will concentrate on this particular case. In scale-free cos- 
mologies two haloes with different masses, M and M' , taken 
at t and t' satisfying the condition a{M,t) = a{M',t') are 
indistinguishable. All the density profiles shown in Figure 2 
correspond to the present time to, but the inside-out growth 
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Figure 4. Same density profiles as in the right panel of Fig. 2 but conveniently shifted so to apply to haloes of the four different masses 
A/; but at four appropriate times t; making them indistinguishable in that scale-free cosmology (provided the assumption of inside-out 
gro wth during accretion). For clar ity we have separated the theoretical profiles derived here (left panel) from the empirical ones obtained 
by Mavarro, Frenk & White (1997) identically shifted (right panel). The shifts applied to the profiles at to can be reconstructed from the 
different symbols used for the non-overlapping empirical profiles: lO^^A/, (full line), 10~^Af, (dot-dashed line), A/« (dashed line), and 
10 Mf (dotted line). Rt{t) is the scale length corresponding to Af, (t). 



condition guarantees that, at any other epoch, haloes have 
the same profiles though truncated at different radii enclos- 
ing the corresponding masses. Thus, by conveniently rescal- 
ing the density profiles of haloes of different masses at to, 
they must perfectly overlap. In Figure 4 we show the result of 
applying such shifts on both the four theoretical density pro- 
files shown in the right panel of Fig ure 2 and the correspond- 
ing empirical profiles obtained by Navarro, Frenk & White 
(19971). I (Shifted density profiles are only displayed down to 
the minimum initial radii where simulated data were fit.) As 
can be seen, the theoretical density profiles satisfy the ex- 
pected homologous condition while the empirical profiles do 
not. There are two possible reasons for this latter fact: either 
the finite resolution of the simulations introduces some ar- 
tificial scale breaking the theoretical homologous behaviour 
of density profiles, or haloes significantly deviate, during ac- 
cretion, from the inside-out growth implied by the adiabatic 
invariant collapse approximation. In the former case, one 
should trust on the theoretical profiles derived here, and 
conversely in the latter. The fact that the deviation of the 
empirical profiles from the homologous condition only con- 
cerns very small radii where the inside-out growth condition 
during accretion should only play an indirect role and where 
the density profiles are the most affected by the finite res- 
olution of the simulations strongly suggests that there may 
be some problem, indeed, with the empirical profiles there. 
It would therefore be of great interest to check through nu- 
merical simulations of scale-free cosmologies whether or not 
large mass haloes grow inside-out during accretion as ex- 
pected from the adiabatic invariant collapse approximation. 

The previous discussion is extremely relevant in con- 
nection with the problem on the actual slopes of density 
profiles at very small radii (r ~ 10~^ R). In general, we find 



good a greement with the slopes reported by Navarro, Frenk 



& Whi ;e (1997) (see Fig. 2). But, for large masses in the 
scale-free case, we find a much steeper slope (as large as —2 



in logarithmic units) . We want to remark that this is neces- 
sary for those density profiles to satisfy the abovementioned 
homologous condition since the slope at small radii of haloes 
of a given mass must match the slope at larger radii of more 
massive haloes. (Note that the asymptotic slope of —2 in the 
scale-free case is independent, indeed, of the value of Am; see 
Fig. 3.) If the correct behaviour of the theoretical profiles de- 
rived here is confirmed we will be led to the conclusion that 
the density profile of very massive haloes deviates from the 
form reported by Navarro and collaborators, and that the 
asymptotic slope at small radii of the halo density profile 
actually depends on the particular cosmology considered. 



3.2 Velocity dispersion profile 

And what about the velocity dispersion profile S(r-)? This 
profile could be readily derived from the Jeans equation (^ 
using the previous density profile provided that the value of 
E were known at some radius. But it is not. (We are avoiding 
of course making any unjustified assumption reporting to in- 
finity.) Since, from the dynamical point of view, the unique 
density profile derived above could perfectly exist side by 
side with an infinite variety of velocity dispersion profiles 
the following questions rises. Why does the empirical veloc- 
ity dispersion profile of haloes show one unique form? What 
fixes it? As mentioned, in the case that the system has un- 
dergone just one unique violent relaxation far in the past and 
no subsequent accretion of matter we would expect it with a 
uniform velocity dispersion. But this solution is incompati- 
ble with the density profile previously calculated imposed by 
accretion. This suggests the following natural guess: the ve- 
locity dispersion profile emerging from violent relaxation at 
the formation of a halo should be as close as possible to uni- 
form, the only small deviation being forced by the boundary 
conditions imposed by accretion at that moment. 

To derive such a guessed velocity dispersion profile we 
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Figure 5. Same as Figure 1 but for the velocity dispersion profile. 







will follow an iterative procedure. For the zero-order solu- 
tion of E'^ at r ^ R{ we take the square value of the uniform 
velocity dispersion that violent relaxation tends to establish, 
given by equation (^) with null logarithmic derivative of E^. 
Then we compute the zero-order value of corresponding 
to a slightly larger or smaller, arbitrarily close radius and 
from the two values we estimate the zero-order logarithmic 
derivative of E^ at r ^ i?f. By substituting this derivative 
into the Jeans equation ^ we obtain the more correct first- 
order value of E^ at that radius. Finally, by repeating this 
procedure, we can obtain an estimate as close as wanted to 
the true value of E^(r) at any radius of the inner region. 
This iterative procedure leads to one unique profile among 
the infinite number that satisfy the Jeans equation for the 
density profile derived above using the mass accretion rate. 
Since we have started the iteration with the values of E(r) 
corresponding to an inner uniform velocity dispersion profile 
(not satisfying that equation) , the final solution is necessar- 
ily the closest to the latter uniform profile which is, at the 
same time, compatible with the boundary conditions im- 
posed by accretion as wanted. 

Once the value of E^(i?f) is known we can use the Jeans 
equation (^) with the universal density profile derived above 
to infer the velocity dispersion profile in the outer region. 
But we can also apply the previous iterative procedure di- 
rectly to any radius of the outer region. Indeed, the fact that 
both solutions satisfy the Jeans equation and that they take 
the same value at R{ guarantees that the two methods yield 
identical results. Moreover, since haloes formed at different 
times with initially relaxed regions overlapping over some 
radial range have, by construction, identical velocity disper- 
sion profiles in the common radial ranges, we are also led to 
the conclusion that there is one unique velocity dispersion 
profile for all haloes having the same mass at a given time 
regardless of their particular formation time, and that this 
unique velocity dispersion profile coincides with the inner 
solution, given by the iterative procedure above, of haloes 
formed at an arbitrarily late time. 

In Figure 5 we compare, for the same cosmologies 
as in Figures 2, the theoretical velocity dispersion pro- 



files, inferred using the previous iterative procedure in the 
whole range of radii, with the e mpirical profiles obtained by 
Navarro, Frenk & White (1997). As can be seen the predic- 



tions are as good as in the case of the density profiles. This 
confirms our guess on the origin of the universal velocity 
dispersion profile. 



4 DISCUSSION 

It is generally believed that the PS formalism describes the 
mass growth of dark-matter haloes in hierarchical cosmolo- 
gies but does not tell us anything about their internal struc- 
ture. This is however in contradistinction with the idea that 
such a structure arises just from gravitation (a scale-free 
force) and the way that dark-matter haloes grow. In the 
present paper, we have shown that, when the distinct dy- 
namical effects of minor and major mergers are taken into 
account, the PS model also makes definite predictions on the 
internal structure of haloes, the resulting average theoretical 
profile being in very good agreement with that drawn from 
high-resolution numerical simulations. Our results therefore 
prove that the structure of haloes, including their scaling 
with mass, is the natural consequence of the combined ac- 
tion of minor and major mergers: outside the total radius of 
the system at formation the structure is fully determined by 
the rate at which mass is accreted, while, inside that radius, 
it results from the initial violent relaxation with boundary 
conditions set by accretion at that moment. The behaviour 
of the theoretical profiles derived here strongly suggests that 
the larger the halo mass, the more apparent are the effects 
on the respective density profiles of the limited resolution 
of the simulations. The confirmation of this suspicion would 
imply a different density profile for very massive haloes from 
that reported by Navarro, Frenk & White (1996, 1997) and, 
what is more important, the non-universality of halo struc- 
ture. 

The aim of the present paper was not to derive the 
accurate distribution of matter in real haloes, but to un- 
derstand the origin of their empirical apparently universal 
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structure. For this reason we have considered the simpli- 
fied case of pure dark-matter haloes as those dealt with in 
A'^-body experiments. In the real universe, about one tenth 
of the halo mass is in the dissipative baryonic component, 
which might have appreciable effects. Likewise, we have as- 
sumed spherical symmetry and neglected any rotation of 
haloes as well as any anisotropy of the velocity tensor while, 
in hierarchical cosmologies, haloes have a slight angular mo- 
mentum and, what is more important, they are immersed 
in large filamentary structures making them a ccrete mat - 
ter preferentially along one privileged direction ( West 1994 ) 
and feel, in their final steady configuration, the tidal field o f 
such anisotropic structures ( Salvador- Sole & Solanes 1993). 



On the other hand, even in the case of exact spherical sym- 
metry, some velocity anisotropy would emerge in the outer 
part of haloes owing to the distinct evolution of the radial 
and tangential velocity dispersions in accreted layers dur- 
ing infall. All these secondary effects explain the elongation 
and slight angular momentum and velocity anisotropy ob- 
served in real as well as simulated haloes. (Note however 
that the density profile derived here is independent of the 
actual degree of anisotropy of the velocity tensor.) Finally, 
in the present paper, we have focused on the structure of 
haloes at z = 0, the only redshift for which accu rate em- 



pirical data for differ ent cosmologies are available (Navarro 



Frenk & : White 1997). In a forthcoming paper, we will study 
the predicted dependence on redshift of this structure and 
compa re it with the results of numerical simulations carried 
out by [Bullock et al. (200l|) for a ACDM cosmology. 
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